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ПРИМЕНЕНИЕ ВЕЙВЛЕТ-ПРЕОБРАЗОВАНИЯ 
И МЕТОДА ПРОНИ ДЛЯ ИДЕНТИФИКАЦИИ БИОГЕННЫХ СИГНАЛОВ 


Рассмотрены нестационарные сигналы биогенного характера (кардиограмма, электроэнцефалограмма 
и др.). Применен анализ сигналов с использованием вейвлет-преобразования и метода Прони, которые 
в отличие от классического метода Фурье дают улучшенную спектральную оценку биогенных сигналов. 
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Введение. Технологии цифровой обработки сигналов и их изображений находят все более ши- 
рокое применение в различных областях, в частности при определении характеристик и иденти- 
фикации биогенных сигналов. 

Разработка алгоритма для анализа биогенного сигнала, однако, является непростой зада- 
чей; довольно часто это даже не целенаправленный процесс. Инженер или компьютерный анали- 
тик часто бывает поражен изменчивостью и разнообразием признаков в биогенных сигналах и 
системах, где эти факторы проявляются в большей степени, чем в физических системах или на- 
блюдениях. Учет всех возможностей и степеней свободы в биогенных системах является наиболее 
сложной проблемой для большинства применений. Методы, показавшие свою пригодность для 
работы с определенными системами или наборами сигналов, могут оказаться несостоятельными 
в других, на первый взгляд похожих, ситуациях. 

Так, для получения правдивого результата биогенного сигнала следует иметь расслаб- 
ленное и спокойное состояние при отсутствии движений. Кашель, напряжение мышц, движения 
глаз, конечностей вызывают соответствующие сигналы, играющие роль нежелательных артефак- 
тов. В результате рекомендацией при снятии сигнала служит задержка дыхания на несколько се- 
кунд. Но это предложение неприемлемо для пациентов в критическом состоянии или при прове- 
дении диагностики у младенцев. 

В таких и других случаях для удаления артефактов используются методы цифровой обра- 
ботки сигналов. 

Обработка сигналов, связанная с анализом их спектров, называется спектральным анали- 
зом. Спектральный анализ используется при распознавании, обнаружении и сжатии сигналов. 
В основе спектрального анализа сигналов лежит интегральное преобразование и ряды Фурье. 

Идентификация особенностей биогенных сигналов — весьма важный этап, так как допу- 
щенные ошибки сказываются на правильности врачебного заключения [1]. 

Постановка задачи. Исследовать методы идентификации для устранения артефактов, способ- 
ных искажать биогенные сигналы, без ухудшения качества исследуемых сигналов. Применить хо- 
рошо разработанные методы Фурье с помощью дискретного преобразования Фурье (ДПФ), ис- 
пользуя алгоритм быстрого преобразования Фурье (БПФ). Применить вейвлет-преобразование и 
метод Прони для анализа и параметризации нестационарных сигналов. Провести оценку погреш- 
ности и сравнительный анализ полученных результатов. 

Методы обработки нестационарных сигналов. Большинство биогенных сигналов имеет 
сложные частотно-временные характеристики. На рис.1 представлен нестационарный кардиосиг- 
нал биогенного характера, на примере которого рассматривается идентификация сигналов. Как 
правило, такие сигналы состоят из близких по времени, короткоживущих высокочастотных и дол- 
говременных, близких по частоте низкочастотных компонент. 
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Рис.1. Кардиосигнал (ЭКГ), содержащий шумы 


Для анализа таких сигналов нужен метод, способный обеспечить хорошее разрешение и 
по частоте, и по времени. Первое условие требуется для локализации низкочастотных состав- 
ляющих, второе — для разрешения компонент высокой частоты. 

Преобразование Фурье - сигнал, заданный во временной области, в виде разложения по 
ортогональным базисным функциям (синусам и косинусам), при этом выделяются частотные ком- 
поненты [2]. Недостаток преобразования Фурье заключается в том, что частотные компоненты не 
могут быть локализованы во времени, что накладывает ограничения на применимость данного 
метода к ряду задач (например, в случае изучения динамики изменения частотных параметров 
сигнала на временном интервале). 

Вейвлеты представляют собой семейства математических функций определенной формы, 
которые локальны во времени и по частоте, и в которых все функции получаются из одной базо- 
вой (порождающей) посредством ее сдвигов и растяжений по оси времени. Вейвлетные функции 
базиса позволяют сконцентрировать внимание на тех или иных локальных особенностях анализи- 
руемых процессов, которые не могут быть выявлены с помощью традиционных преобразований 
Фурье. 

Метод Прони - это метод анализа коротких отрезков сигнала, основанный на аппроксима- 
ции сигнала конечной суммой комплексных экспонент (рис. 2-3), т.е. на подгонке экспоненциаль- 
ной модели к измеренным эквидистантным (равноотстоящим, равноудаленным) значениям и по- 
следующем вычислении дополнительных значений посредством оценки параметров этой экспо- 
ненциальной модели в промежуточных точках. 
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Рис. 2. Короткий отрезок кардиосигнала, 
обработанный методом Прони (сумма комплексных экспонент) 
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Рис. 3. Спектр короткого отрезка сигнала 
и фаза после анализа методом Прони 


Краткий обзор преобразования Фурье. Классическим методом частотного анализа сигналов 
является прямое преобразование Фурье: 


©. © 
5(@)= [5(е "4. (1) 
— 00 
Обратное преобразование Фурье: 
о = | $(®)е/^'а®. (2) 
21 


Результат преобразования Фурье — амплитудно-частотный спектр, по которому можно оп- 
ределить присутствие некоторой частоты в исследуемом сигнале (рис.4). 
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Рис. 4. Фурье-спектр сигнала 


Преобразование Фурье имеет ряд недостатков: 
— исходный сигнал заменяется на периодический с периодом, равным длительности ис- 


следуемого сигнала; 
— со временем плохо работает при изменении параметров процесса (нестационарности), 


поскольку дает усредненные коэффициенты для всего исследуемого образца; 
— не дает представления о динамике изменения спектрального состава сигнала. 
Рассмотрим три подхода к анализу нестационарных сигналов биогенного происхождения — 


метод оконного преобразования Фурье, метод вейвлет-преобразования и метод Прони. 
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Оконное преобразования Фурье. Классическое преобразование Фурье связано со спектром 
сигнала, взятым во всем диапазоне существования переменной. При анализе нестационарных 
сигналов интерес представляет только локальное распределение частот, в то время как требуется 
сохранить изначальную переменную (обычно время). В этом случае используется оконное преоб- 
разование Фурье, представляющее собой сегментирование сигнала на фрагменты («окна»). Для 
начала необходимо выбрать некоторую оконную функцию: 


Е( о) = р | хфи(-де““ат, (3) 


где Е( ©) — распределение частот (несколько искаженное) части оригинального сигнала /(Г) 


в окрестности времени {. 
Недостатки оконного преобразования Фурье. Когда нестационарный сигнал имеет высоко- 
частотные компоненты в течение короткого промежутка времени и низкочастотные колебания 
при рассмотрении больших временных масштабов, оконные преобразования позволяют проана- 
лизировать либо высокие частоты в коротком окне времени, либо низкочастотную компоненту, но 
не оба колебания одновременно (рис.5). 
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Рис.5. Кардиосигнал после обработки методом Фурье 


В результате был разработан подход, при котором для различных диапазонов частот ис- 
пользовались временные окна различной длительности. Оконные функции получались после рас- 
тяжения-сжатия и смещения по времени гауссиана. Эти базисные вейвлет-функции называются 
компактными волнами. 

Функция ту (1) может называться вейвлетом, если она удовлетворяет следующим свойствам: 

— область определения функции - вся числовая ось; 

— интеграл от функции ту (#) по области определения равен нулю; 

— функция чу (1) быстро убывает при ЕЁ >-00. 

На рис.6 изображены примеры простейших базисных функций, посредством которых могут 
быть представлены произвольные колебания. Вейвлет-разложение по ним проводится путем вы- 
числения сверток сигнала /` с компактной волной ‘у при различных масштабах и сдвигах аргумента. 
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Рис.6. Компактные волны, полученные дифференцированием функции Гаусса 
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Рассмотрим свертку произвольного сигнала / (1) с функцией: 
и (та) = Иа | уу (1— </а 1, (4) 


где т- дрейфовая переменная; а — масштабная переменная. 

Результатом свертки является функция, определенная в частотно-временной области, 
в отличие от преобразования Фурье, которое определено только в частотной области. Таким об- 
разом, вейвлет-преобразование позволяет судить не только о частотном спектре сигнала, но так- 
же о том, в какой момент времени появилась та или иная гармоника. 

Параметр //а — частотная характеристика сигнала, подвергшегося вейвлет-преобразова- 
нию, привязана к «временной точке» - т. Это означает, что с помощью вейвлет-анализа мы мо- 
жем изучать сигналы, частотные характеристики которых локально зависят от времени. Фурье- 
анализ такой возможности не дает. 

Как и Фурье, вейвлет-преобразование имеет обратное преобразование, которое записы- 
вается формулой: 


00 со 


Х(х)=з ус -Иа) (Ба) ааа? СЕ, (5) 


—о 0 
где © — некоторая константа [5]. 
Вейвлет-преобразование Добеши (рис.7). 
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Рис.7. Структура вейвлета Добеши 
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С _С а С — матрица Добеши. 

С С | | | о - : с с Преобразование Добеши от неко- 
с С И С ь торого вектора соответствует ум- 
о о ножению его на матрицу. 


Каждая строка матрицы соответствует свертке вектора сигнала с локализованным вейв- 
лет-фильтром. При этом нечетные строки сглаживают, а соответствующие четные фильтры вы- 
числяют поправки к сглаживанию, обеспечивая обратимость преобразования. Коэффициенты 


С, ...С. могут быть выбраны так, чтобы поправки обращались в нуль для гладких функций. За- 


459 


Физико-математические науки 


нуление свертки с постоянной и линейной функциями дает два условия на неизвестные коэффи- 
циенты: 
= =, =0 
ОС; —-1С> +2С; -3С =0. 
Далее потребуем, чтобы преобразование было ортогональным — обратное преобразование 


получалось простым транспонированием. Из условия равенства произведения матриц прямого и 
обратного преобразования единичной матрице получаем еще два условия [4]: 


С? + С? + С? + С? =1 >  @=(1+\3)/4]2 с, = (3-\3)/4\/2 
СС +СС =0 Откуда с, = (3+\3)/42 су = (1-3) /4\/2. 


Полное вейвлет-преобразование получается путем умножения вектора сигнала на полу- 
ченную матрицу Добеши, далее нечетные компоненты сигнала (сглаженная версия сигнала) снова 
умножаются на матрицу Добеши (вдвое меньшей размерности), и так до тех пор, пока не будет 
получено полностью сглаженное представление, т.е. среднее значение сигнала (рис.8-9). 
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Рис.8. Кардиосигнал, обработанный вейвлет-методом 
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Рис.9. График кардиосигнала и спектрограмма 


Исходный подход Прони. Если число используемых отсчетов данных равно числу экспоненци- 
альных параметров, то возможна точная подгонка экспонент под имеющиеся данные. Рассмотрим 
функцию дискретного времени, представляющую собой сумму р-экспонент: 
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Р 
Жи= У Ир. (6) 
= 


Заметим, что в этом выражении используется х[и], а не х [п], поскольку точно 2р-ком- 
плексных отсчетов х[1], ..., х[2р] используется для точной подгонки к экспоненциальной модели 
с 2р-комплексными параметрами Й,, ...,Й», 2/1 ..., 2». Входящие в р-уравнения (6) можно записать 
в матричной форме: 


0_0 0 
2] 2) ИЕ 2р й ХИ] 
212) све 2. й> м2] 
= || (7) 
р В й, МР] 


Матрица с временными индексами элементов 2 имеет структуру матрицы Вандермонда. 
Если может быть найден метод для раздельного определения элементов 2, то уравнение (7) мож- 
но рассматривать как систему уравнений, решив которую, определяют неизвестный вектор ком- 
плексных амплитуд. 

Ключ к разделению основан на том факте, что уравнение (7) является решением некото- 
рого однородного линейного разностного уравнения с постоянными коэффициентами. Для того 
чтобы определить вид этого разностного уравнения, определим сначала полином ф(2), корнями 


которого являются экспоненты 2х. 
Р 
(2) =] [ (2-2%). (8) 
К=1 


Если произведения в (13) выразить в виде степенной последовательности, то полином 
можно представить в следующем виде: 


р 
(=) =>. атреР", (9) 
т=0 
с комплексными коэффициентами а[т], для которых а[0]=|. 
Осуществляя в уравнении (6) сдвиг индекса от и к (и-т) и умножая обе его части на па- 
раметр а[т], получаем 


р 
а[т]х[и — т] = ат] У Икар". (10) 
&=1 
Записывая аналогичные произведения а[О]х|[и], ...‚а[т-Пх| и-т-+1] и осуществляя сумми- 
рование, получаем: 


р р р 
уз ат]х[п — т] = >. й. рэ ат]2" "1 | (11) 
т=0 1=0  т=0 

которое справедливо при р+1 <и <2р. Получаем: 

р р р 
У ати —т] = ре" У’атР" =0. (12) 
т=0 1=0 т=0 


Сумму в правой части (12) можно рассматривать как полином, определяемый уравнением 
(9), который записан через свои корни, что и обеспечивает в уравнении (12) равенство нулю. 
Уравнение (12) - это линейное разностное уравнение, однородное решение которого выражается 
формулой (6). Полином (9), ассоциированный с этим линейным разностным уравнением, называ- 
ется характеристическим. 
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Можно записать в виде следующего рхр-матричного уравнения р-уравнения, представ- 
ляющие истинные значения коэффициентов а[и|, удовлетворяющих (12): 


ЖРХР-И........щ... хи у а хр-+П 
МР + ИР] елььиь, х2] | а[2] мр-+2] 

=. (13) 
хр - Пар -2]..Мр ар х[2 р] 


Из уравнения (13) следует, что, имея 2р-отсчетов комплексных данных, возможно разде- 
ление множеств параметров й/ и 2„. Комплексные полиномиальные коэффициенты а[1], ..., а[р], 
которые являются функциями только зависящих от времени компонентов 2, экспоненциальной 
модели, позволяют по временным отсчетам сформировать соотношения для линейного предска- 
зания. Матрица в уравнении (13) имеет тёплицеву структуру [3]. 

Процедуру Прони для подгонки р-экспонент к 2р-отсчетам данных можно теперь предста- 
вить в виде следующих трех этапов. 

На первом этапе получается решение уравнения (13) для коэффициентов полинома. 

На втором этапе вычисляются корни полинома, определяемого уравнением (9). Используя 
корень 2; можно определить коэффициент затухания а; и частоту синусоиды /; с помощью соотно- 
шений: 


а = ш|2;|/Т(с'), (14) 
Х. = агсю [п {2, /Ве{2, Н/2^Т(Ги). (15) 
Для завершения процедуры Прони корни полинома, вычисленные на втором этапе, ис- 
пользуются далее для формирования элементов матрицы уравнения (7), которое затем решает- 
ся относительно р комплексных параметров 1[1], ...., #[р]. Каждый параметр 7; используется да- 
лее для определения амплитуды 4; и начальной фазы 0;, которые вычисляются с помощью ва- 
риаций 

0, = агсёе [ша {й; }/ Вей, й (рад). (17) 

Результаты вычислений показаны на рис.3, 10-11. 
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Рис.10. Кардиосигнал, обработанный методом Прони 
(сумма комплексных экспонент) 
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х 10 АЧХ-метод Прони 
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Рис.11. Спектр сигнала и фаза после анализа методом Прони 


Метод наименьших квадратов Прони. На практике число отсчетов данных Л, как правило, пре- 
вышает то минимальное их количество, которое необходимо для подгонки модели из р-экспонент, 
т.е. № > 2р. В этом переопределенном случае последовательность отсчетов данных может быть 
аппроксимирована лишь как экспоненциальная последовательность: 


р 
х[п] = У ик | ‚ ге |1 <п <^№. (18) 
К=1 
Ошибка аппроксимации в данном случае распределяется выражением =[и] = х[п]- х[п]. 


Одновременное нахождение порядка р и параметров {й,, 2„}, где 1 < К < р, которые минимизируют 
сумму квадратов ошибки и представляет собой трудную нелинейную задачу: 


М 
р=У Ня]. (19) 
п=1 


Используя вариант метода Прони можно определить субоптимальное решение, которое 
обеспечивает удовлетворительные результаты. Используя на первом и втором этапе трехэтапного 
метода Прони соответствующие линейные процедуры наименьших квадратов, получим процедуру 
экспоненциального моделирования, которую иногда называют обобщенным методом Прони [3]. 
Выводы. Применены хорошо разработанные методы Фурье с помощью дискретного преобразо- 
вания Фурье (ДПФ), использовался алгоритм быстрого преобразования Фурье (БПФ), вейвлет- 
преобразование и метод Прони для анализа и параметризации нестационарных сигналов. Прове- 
дена оценка погрешности и сравнительный анализ полученных результатов. 

На основе результатов анализа ЭКГ человека показано, что вейвлетные базисы могут 
быть хорошо локализованными как по частоте, так и по времени. При выделении в сигналах хо- 
рошо локализованных разномасштабных процессов можно рассматривать только те масштабные 
уровни разложения, которые представляют интерес. 

Вейвлетные базисы, в отличие от преобразования Фурье, имеют достаточно много разно- 
образных базовых функций, свойства которых ориентированы на решение различных задач. Ба- 
зисные вейвлеты могут иметь и конечные, и бесконечные носители, реализуемые функциями раз- 
личной гладкости. 

Использование метода Прони позволяет получить более точные спектральные оценки для 
коротких последовательностей данных. Указанные возможности позволяют получить новую ин- 
формацию о динамике изменений ЭКГ в процессе деятельности человека. 
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